Light scattering for thermometry of fermionic atoms in an optical lattice 
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We propose a method for measuring the temperature of fermionic atoms in an optical lattice 
potential from the intensity of the scattered light in the far-field diffraction pattern. We consider a 
single-component gas in a tightly-confined two-dimensional lattice, illuminated by far off-resonant 
light driving a cycling transition. Our calculations show that thermal correlations of the fermionic 
atoms generate fluctuations in the intensity of the diffraction pattern of light scattered from the 
atomic lattice array and that this signal can be accurately detected above the shot noise using a 
lens to collect photons scattered in a forward direction (with the diffraction maxima blocked). The 
sensitivity of the thermometer is enhanced by an additional harmonic trapping potential. 

PACS numbers: 03.75.Ss,42.50.Ct 



Ultra-cold atomic gases in optical lattices can consti- 
tute almost ideal realizations of Hubbard models 0, Q 
that are fundamental to strongly-correlated physics. Re- 
cent experiments on fermionic atoms in lattices have 
demonstrated both a superfluid pairing Q and a Mott 
insulator 0, Q, opening up possibilities for experimen- 
tal simulation of even more complex strongly-correlated 
systems, such as antiferromagnetic phases and high-T c 
superconductivity. Thermal energy is a major control 
parameter in fermionic lattice systems and characteriz- 
ing phase diagrams is fundamentally related to the ability 
to perform accurate temperature measurements. 

Here we show that the temperature of fermionic atoms 
in a 2D optical lattice can be accurately determined by 
measurement of the light scattering from the atoms. The 
diffraction pattern is insensitive to thermal atomic corre- 
lations but blocking the diffraction maxima and collect- 
ing light scattered outside the diffraction orders using 
a lens provides a measurable optical signal that reflects 
thermal and quantum fluctuations of lattice atoms. 

The temperature of fermionic atoms in optical-lattice 
experiments has been deduced indirectly from the tem- 
perature of the trapped cloud of atoms before turning up 
the lattice 0, Q , and by detecting double-occupancy in a 
two-species gas by converting atom pairs into molecules 
[3, 0| ■ Other temperature measurements detected atomic 
shot-noise Q or the sharpness of interference peaks Q in 
absorption images after a ballistic expansion. This exist- 
ing technology provided vital information about temper- 
ature but it has limitations, and there is a clear need for 
new methods; e.g., detecting atomic shot noise proved 
inconclusive in some superfluid/thermal lattice systems 
Q and it was argued that detecting superfluidity from 
the sharpness of interference peaks can be ambigious as 
even a thermal gas may show misleadingly sharp peaks 
Q. Moreover, a range of phenomena occur below the 
Neel temperatures [l0( of these systems, e.g., antiferro- 
magnetic ordering and superfluid pair hopping, but this 
requires significantly more cooling than current exper- 
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FIG. 1: (Color online) Left: The arrangement for diffraction 
from a 2D optical lattice. The incident light propagates in 
the positive z direction and scatters off a regular array of 
atoms in the xy plane. The two lenses have focal lengths 
/i and /a. Right: View along the z-direction showing the 
intensity pattern of elastic scattering in the focal plane of lens 
1 (with log scaling, as in Fig. [2}, and the square (A), or cross- 
shaped stop (B), used to block the incident beam and most 
of elastically scattered light in central (zero order) diffraction 
peak (higher orders fall outside the range of angles collected 
by the lens). In the focal plane of lens 2 a photodiode, or 
similar detector, measures the intensity of the (unblocked) 
scattered light. 



iments [4j further increasing the need for a method to 
measure ultra-low temperatures. The diffraction of light 
from regular arrays of atoms in an optical lattice re- 
sembles the powerful method of x-ray diffraction from 
crystalline materials. The diffraction pattern reflects the 
atomic lattice structure and the overall diffraction en- 
velope the shape of atomic wavefunctions at individual 
sites. Scattering into angles outside the diffraction orders 
arises by inelastic processes in which the phonon excita- 
tions of atoms in the lattice absorb the recoil kicks from 
the scattered photons and generate fluctuating shifts in 
the diffraction pattern which carry information about 
thermal and quantum correlations of the atoms. Light 
scattering from optical lattice systems was previously 
studied theoretically for detecting particle number frac- 
tionalization (ll| . coupling of atoms to optical cavities 
[l2| and to polarization of light [l3| ■ 

We consider a fermionic atomic gas in an optical lattice 
illuminated by light with the positive frequency compo- 
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ncnt of the electric field amplitude E^(r) = i£e; n e lk ' r , 
the polarization e; n , and the wave vector k = fce z ; see 
Fig. [TJ In lattices the atom dynamics can be restricted 
to ID (highly-elongated tubes) or 2D (pancake-shaped 
layers) by optical confinement For light passing 

through in an appropriate direction such a sample is opti- 
cally thin and this makes light scattering a more suitable 
probe than in 3D samples. In the limit of a large fre- 
quency detuning 6 of the light from the atomic resonance 
(compared to the Weisskopf-Wigner linewidth 7), the dy- 
namics of the electronically excited atomic state may be 
adiabatically eliminated [15j. The scattered light ampli- 
tude E+(r) then becomes proportional to the transition 
amplitude of atoms between the initial and final hyper- 
fine (electronic ground) states, g and g', and we calculate 
the scattered light at r (in the far radiation field lion by 
integrating over the atomic dipoles residing at r' jl5| |: 



E+(r) = CA g , g 



d 3 r / e -iAk.r', 



*;,(r')* fl (r'). (1) 



Here Ak = k' — k ~ k(h— e z ) is the change of wave vector 
of light upon scattering (|Ak| = 2k\ sin(0/2)|) with n be- 
ing a unit vector in the direction of the detector from the 
sample, C = 3^e lkr j /(ASkr). A g > g = n x (n x d s < e )(e in • 
d eff )/D 2 , and d eg = (d ffe )* = (e|d|.g) = D(e\ag)e* a is the 
atomic transition dipole matrix element between g and 
an excited Zeeman state e, where D = (Gnheo^/k 3 ) 1 / 2 
is the reduced dipole matrix element and (e\ag) are the 
Clebsch-Gordan coefficients. Here we use the convention 
that repeated indices (3, e,tr, etc.) are summed over. 

Ultra-cold atoms only occupy the lowest vibrational 
state of each lattice site and we can expand the mat- 
ter field in Eq. ([1]) in terms of the Wannier functions 
w g j (r) = w g (r — Tj ) localized at each site j centered at 
Tj, i.e., we write ^(r) = J^j w g,j( r )bgj- Here b g j denotes 
the corresponding annihilation operator. We obtain for 
the scattered light intensity I(r) = 2e c(E sc E+) which 
carries information about the atomic correlations 
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Here M^w 



ii„(37/2<5fcr) 2 with 
^in = £oc£ 2 /2 being the incoming light intensity. The 
Debye- Waller factor a(h) depends on the integrals of w g j 
and, for the case of a level-independent lattice potential, 
it is obtained from the Fourier transform of the lattice 
site density a — \ J d 3 r e~ iAk ' r |w s (r)| 2 | 2 . 

Equation ([2j) gives the scattered light intensity for an 
arbitrary lattice system if the appropriate atomic corre- 
lation functions can be evaluated. Here we shall apply 
it to a single-component, non-interacting FD gas in a 
2D square lattice, as illustrated in Fig. Q] — inelastic scat- 
tering events that carry information about thermal and 
quantum fluctuations of atoms are mapped onto the light 
field, hence providing a sensitive in situ thermometer. 



We assume that atoms initially occupy only one hy- 
pcrfine level and that atoms scatter back to the same 
level (cycling transition). Consequently, we drop any 
explicit reference to hyperfinc levels g. We consider 
a cr + -polarizcd light beam, with wavelength 766.5 nm, 
driving |1) = \AS 1/2 ;F = 9/2, m F = 9/2) -> |2) = 
|4P 3/2 ; 11/2, 11/2) transition of 40 K, so that d 12 = De +1 
and M\'\ = (3 + cos2#)/4. In a single-component 
fermionic gas s-wave scattering between ground state 
atoms is forbidden and, at low temperatures, the atoms 
can be considered as non-interacting. 

We consider a uniform 2D square lattice with periodic- 
ity a, potential V = sEr [sin 2 (7ra;/a) + sin 2 [iry / a)] , and 
the lattice-photon recoil energy Er = ir 2 H 2 / (2ma 2 ). The 
atom dynamics along the z-axis is assumed to be frozen 
out by strong confinement. The Wannier functions Wj(r) 
can be approximated by ground state harmonic oscillator 
wave functions with frequencies ui z and (j XtV = 2^fsEf>Jh, 
obtained by expanding the potential at the lattice site 
minimum. Then a(n) = U.i= x , v ,z exp{ — ((Afc i ) 2 ^/2)}, 
with I4 = (h/muji) 1 / 2 . The only contribution to the 
Hamiltonian H = — J^iplbj +bjbi) arises from the hop- 
ping J ~ 4s 3 / 4 e _2v ^i?fl/v / 7F, where the summation is 
over adjacent sites only. 

The Hamiltonian is diagonalized in the quasi- 
momentum space by bj = (1/N S ) ^ q u q a q e lq rj ', where 
the amplitudes |u q | = 1, and N s is the number of 
sites in each direction (both x and y). The energy 
is given E^ = 4J {sin 2 (q x ja/2) + sin 2 (q y jCi/2)}, with 
qj = (qx 3 ,qyj) = (2ir/N s a)(j x ,j y ), where the integers 
may be chosen as j x<y = —N s /2, . . . , (N s /2 — 1). We 
may now calculate the intensity in Eq. ^ by evaluating 
(b\biP-bj) using the phonon operators, a q , to give 



I/A = / 2 2l A - k + E ~ 21 



Ak+q'-q ! 
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where A = M{l\a(n)B, f = N/N 2 < 1 is the filling 
fraction of the lattice with TV being the total number of 
atoms, n q denotes the occupation numbers of the ideal 
FD distribution, and 
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(N s Ak ja /2) 
(Afcja/2) 



(4) 



is the diffraction pattern from a 2D square array of _/V s x 
N s diffracting apertures; where Ak is the change in the 
wave vector of light on the xy-plane. 

The first term in Eq. ©, proportional to f 2 , is the elas- 
tic scattering contribution where an atom scatters back 
to its original cm. state. It generates the diffraction pat- 
tern from a non-fluctuating atom density; see Fig. [TJ For 
a large number of sites, this gives rise to sharp diffraction 
peaks of small angular spread with significantly weaker 
intensity between the orders. The second term in Eq. ([3]) 



3 




FIG. 2: (Color online) The angular distribution of the in- 
tensity (in units of B defined under Eq. ((2J)) for elastic 
scattering (left on a logarithmic scale) at scattering angles 
4> — 0, 0.025, 0.1 rad (top to bottom); and for inelastic scat- 
tering (right on a linear scale) showing the effect of fermionic 
statistics at temperatures T/T F = 0.01,0.05,0.25,0.5 (bot- 
tom to top) with Tf = 64 nK and = 0. Elastic scattering 
falls off faster in directions away from the x and y axes (Fig.[T]) 
therefore quantitative comparison with the inelastic scatter- 
ing requires integration over the collected angles (Fig. [3]). 
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FIG. 3: (Color online) Left: The total number of collected 
photons/second (elastic + inelastic) as a function of temper- 
ature for a square stop (solid line), a cross-shaped stop of 
angular width 0.01 rad (dashed line) and 0.05 rad (dashed- 
dotted). Right: The required number of experimental rep- 
etitions needed to measure temperature to certain relative 
precision (AT/T2), with solid (dashed) lines for the square 
(0.01 rad cross-shaped) stop and temperatures of T/Tf = 
0.02, 0.05, 0.15, 0.25 (top to bottom). 



represents inelastic scattering where an atom is excited 
from a quasi-momentum state q and scatters to a differ- 
ent state q'. It is this term that contains the effect of 
the FD statistics, in both (a) the quasi-momentum dis- 
tribution fiq that obeys the quantum statistics, and (b) 
the product of the occupation numbers that describes 
the Fermi inhibition: scattering events in which an atom 
would recoil into an already occupied state are forbid- 
den by the Pauli exclusion principle. Importantly, this 
inelastic term leads to scattering outside the diffraction 
orders — phonon excitations lead to a deflection of the 
light since an atomic recoil Aq = q' q is associated with 
an equal and opposite change of photon momentum. In 
a sequence of experimental realizations the values of n q 
fluctuate according to the FD statistics. Therefore the 
second term in Eq. ([3]) generates, on the top of the diffrac- 
tion pattern, an inelastically scattered light intensity that 
also fluctuates during the sequence measurements. 

To illustrate this method, we study light scattering 
from a N s = 150 lattice of 40 K atoms with / = 0.5, 
s = 7.8, J = 0.04 E R , a = 0.4 /mi, and </> = tt/2. The re- 
coil component of the photon on the xy-plane is absorbed 
by the atom quasi-momentum with | Aq| = k sin((9). Pho- 
ton recoil to higher bands on the xy-plane and in the 
z-direction is negligible (because we have taken hu z ~ 
7.5E' R , where E' R is the probe beam recoil energy, cor- 
responding to l z ~ 63 nm and u z /2tt ~ 63 kHz, and the 
energy band gap on the zy-plane is about 5E' R when 
the maximum energy absorbed by atoms due to a recoil 
kick in the xy-plane is E' R (l7]). At T = the atoms 
fill the Fermi sea n q = Q(qF — |q|) and only inelas- 
tic scattering events for which the final state is out of 
the Fermi sea are allowed. Thus all scattering events for 
which it /a > |Aq| > 2qp are allowed, but for small an- 
gles sin(0) < 2qp/k some recoil events would lead to an 
already occupied state and are forbidden [lij . Thus in- 
elastic scattering is strongly suppressed by the FD statis- 
tics in the near-forward direction (or near the diffraction 
maxima); here 2qp/k ~ 1.36 > 1 and it is at least par- 



tially suppressed in all directions as seen at low T in Fig. [2] 
(right). At large angles (9 > 0.5 rad) the diffraction pat- 
tern envelope is determined by the internal atomic level 
structure and the lattice site wave function (hence does 
not give information about the FD statistics). Elastic 
scattering exhibits a diffraction peak in the forward di- 
rection which we block using a suitable stop (Fig. [T|), 
together with any unscattered part of the incident beam; 
this light does not provide any information about the 
atom statistics. 

In Fig. [3] we show the number of collected pho- 
tons/second by the lens with the numerical aperture 
NA=sin0 max =O.48 for I m = 5Wm~ 2 and S = 2O7. 
We use a square stop of half-diagonal angular width of 
0.05 rad with its diagonals oriented along the x and y 
axes, and a cross-shaped stop, \6 X \ < 0, \9 V \ < <d, for 
G = 0.01 and 0.05 rad to block most of the elastically 
scattered light without dramatically reducing the detec- 
tion of inelastically scattered photons [3] . 

In each inelastic scattering process an atom absorbs a 
recoil kick which changes its cm. state. This affects the 
temperature (heats the sample) and limits the maximum 
number of inelastic scattering events W in a single ex- 
perimental realization of the lattice system to be a small 
fraction W/N of the total number of atoms (otherwise the 
measurement perturbs the ground state significantly). 
Only a fraction rj(T) of all the inelastically scattered pho- 
tons are collected by the lens, so N™(T) = -q(T)W is 
the number of detected inelastically scattered photons 
in each realization (assuming 100% detector efficiency). 
In experiments the lattice system can be prepared and 
measured r times so that the total number of detected 
photons is tN c (T) = r[iV,j n (r) + iV c cl (T)], where N? de- 
notes the number of detected, elastically scattered pho- 
tons in a single realization. For far off-resonant light, the 
fluctuations of the number of scattered photons are Pois- 
sonian, so to distinguish between two different optical re- 
sponses corresponding to temperatures If and Ti (with 
T2 > T\), the difference between the number of detected 
photons in the two cases must be of order \J tN c (T2 ) . 
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FIG. 4: (Color online) The angular distribution of the scat- 
tered light intensity (at = 0) from a ID lattice (left), 
and the same lattice with an additional harmonic poten- 
tial in which lh = 2 fim (right). The rapidly oscillatory 
curves on both sides represent elastic scattering (on right 
for T = 0.5 TV) and the smooth curves correspond to in- 
elastic light at T/T F = 0.01,0.05,0.25,0.5 (bottom to top at 
9 = 0). The inset shows the temperature dependence of the 
number of collected photons/second (elastic + inelastic) for 
a stop with \9 X \ < Q for the two cases: no trap, NA=0.48, 
O = 0.047 rad (dashed line); and atoms in a potential with 
l h = 2/xm, NA= 0.8, 9 = 0.086 rad (solid line). 



case (N s = 150 with the lattice along the x axis and 
ly,z = 63 nm) we find that the angular distribution of 
the scattered light for a trap with lh = 6 /J,m is modi- 
fied very slightly (less than 2% difference at any point), 
whereas the results for an additional potential in which 
lh = 2 /mi are shown in Fig.[¥] The temperature variation 
in the trap is notably stronger than for a translationally 
invariant lattice because the inelastic scattering depends 
on temperature at all angles. Thus the temperature sen- 
sitivity of the measured signal is enhanced. 

In conclusion, wc propose an efficient optical ther- 
mometer for fermionic atoms in a lattice. The general 
scattering formula ((2|) may easily be applied also to multi- 
species systems @, provided that the atomic corre- 
lation functions can be evaluated. 

ABD (via QIP IRC) and JR acknowledge financial sup- 
port from the EPSRC. 



Thus the number of repetitions of the measurement re- 
quired for an uncertainty of AT = T 2 — T\ in the deter- 
mination of temperature is the value of r which satisfies 
t{N c (T 2 ) - AT c (Ti)) ~ y/TN c (T 2 ). 

For simplicity, we set W/N = 0.1. Figure [3] shows 
the number of experimental realizations of the lattice 
system required to achieve an accuracy AT/T 2 , e.g., 
T/Tp = 0.1,0.25 can be measured within 5% by prepar- 
ing the system 98 and 14 times respectively, providinga 
highly sensitive thermometer for atoms [20j. In Ref. [4j 
a Mott insulator state was observed at T ~ 0.28Tp. The 
fraction r){T) of inelastically scattered photons collected 
by the lens is large in Fig. [3] (varying from 18 to 37%). 
This arises because (a) the intensity in the forward di- 
rection is twice the intensity in the transverse direction 
for cr-polarizcd light (a well-known feature of the Zecman 
effect); and (b) the size of the lattice site wave function 
generates an overall envelope for the diffraction pattern 
that suppresses radiation at large angles. We found that 
the signal-to-noise ratio is not improved with a larger 
aperture lens (NA=0.75), but would be increased by us- 
ing a stack of several identical 2D lattice layers, provided 
that the system remains optically thin. 

In typical experiments atoms are confined in an opti- 
cal lattice plus harmonic trap — to check that this does 
not render our method inaccurate we carried out cal- 
culations for a ID lattice system with a harmonic po- 
tential of frequency il. This affects the plane wave ba- 
sis used in Eq. ^ when the variation of the trap en- 
ergy over the sample is greater than the hopping en- 
ergy: mVl 2 (N s a) 2 > J; which for our parameters oc- 
curs when f2 > 27r x 7 Hz (lh = ^h/mQ < 6 fxva) . In a 
trapping potential we diagonalizc the Hamiltonian H = 

E 3 -[Cj 2 && " Jtybj+i +f J+ ih)], with C = (a 2 /^) 2 £fl, 
and use the eigenfunctions as a new basis to evalu- 
ate Eq. ([2]). For the same parameters as in the 2D 
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